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The latest Planck results reconfirm the existence of a slight but chronic tension between the 
best-fit Cosmic Microwave Background (CMB) and low-redshift observables: power seems to be 
consistently lacking in the late universe across a range of observables (e.g. weak lensing, cluster 
counts). We propose a two-parameter model for dark energy where the dark energy is sufficiently 
like dark matter at large scales to keep the CMB unchanged but where it does not cluster at 
small scales, preventing concordance collapse and erasing power. We thus exploit the generic scale- 
dependence of dark energy instead of the more usual time-dependence to address the tension in 
the data. The combination of CMB, distance and weak lensing data somewhat prefer our model to 
ACDM, at = 2.4. Moreover, this improved solution has erg = 0.79 ±0.02, consistent with the 
value implied by cluster counts. 


I. INTRODUCTION 

The recent Planck cosmology results [1, 2] provide 
stunning support for the ACDM “standard model” of cos¬ 
mology. One of the few results that are not in quite as 
excellent agreement with the parameter constraints from 
the measurements of the anisotropies in the cosmic mi¬ 
crowave background (CMB) is the determination of the 
amplitude of cosmological perturbations at late times on 
small scales. This slight tension is most apparent when 
comparing to weak lensing measurements as provided by 
CFHTLenS, which prefers lower values of cts (Ref. [2, 
Fig. 18]) and also the lower-than-expected cluster abun¬ 
dances (Ref. [3, Fig. 10] and [4]). In addition, dark mat¬ 
ter growth rates, proportional to erg) tend to fall on the 
low side of the concordance values, even if they are not 
inconsistent [2, Fig. 16]. Although it is well possible that 
these discrepancies come from systematic effects in the 
analysis of the data in the low-redshift universe, it is also 
the case that late-universe measurements of erg consis¬ 
tently imply a smaller value than that obtained by pro¬ 
cessing the initial amplitude of perturbations observed in 
the CMB through concordance gravitational collapse. 

One possible resolution of this conflict is to suppress 
the clustering at low redshift through a new physical 
mechanism active only at late times. For example, heavy 
neutrinos could play such a role for clusters [5, 6]. How¬ 
ever, the latest data do not seem to support this [2]. 

Somewhat more exotically, dark energy could be dy¬ 
namical, changing the evolution of the universe at late 
times. However, the analysis in [7, Fig. 4] shows that 
distance data from baryon acoustic oscillations (BAO) 
constrains the usual quintessencelike models too much 
to significantly improve the possible fits to weak lens¬ 
ing data. More general models of modified gravity in 
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the quasistatic limit that do not contain ghosts typically 
serve to increase growth rates [8-15], although a tempo¬ 
rary suppression around redshift 0 « 1 is also typical [16] 
and could possibly be exploited with some tuning. 

There is another possibility: dynamical dark energy 
could affect dark-matter clustering in a scale-dependent 
manner. The CMB is mostly a large-scale observable, 
while galaxy weak lensing has only been measured on 
small scales; clusters are also small-scale phenomena. A 
model where DE causes the dark matter to cluster slower 
at small scales than at large without significantly mod¬ 
ifying the background expansion history could produce 
the desired phenomenology. 

We propose a mechanism to achieve this: we exploit 
the dark degeneracy [17] to trade some of the dark matter 
for dark energy. Provided the two are similar enough, ob¬ 
servables remain unchanged with respect to ACDM. We 
investigate a very minimal extension of the concordance 
cosmological model: we allow the dark energy to be a per¬ 
fect fluid with constant pressure and a constant nonzero 
sound speed. This keeps the expansion history exactly 
the same as ACDM. Such a choice allows us to concen¬ 
trate on investigating the effect of changing the proper¬ 
ties of the perturbations of dark energy while keeping 
the background, which is already very well constrained, 
fixed. 

If the dark energy has exactly zero sound speed, then 
the dark degeneracy also ensures that the linear per¬ 
turbations in the metric are unchanged with respect to 
ACDM. Effectively, the dark energy component that re¬ 
places part of the dark matter also clusters like that miss¬ 
ing dark matter. However, if the sound speed is nonzero 
then the dark energy clusters less than the missing dark 
matter on scales inside its sound horizon. We show that 
this model, as a result of this change in the behavior of 
DE at its sound horizon, indeed goes some way to relieve 
the tension between the early and late universe and is 
favored over ACDM when the combination of CMB and 
weak lensing data are used. Moreover, the best-fit model 
has a lower value of as which should help with the tension 
with cluster number counts. 
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II. DARK-ENERGY MODEL 

All measurements of the cosmological background ex¬ 
pansion history, such as baryon acoustic oscillations 
(BAO) or supernovae (SNe), or even cosmic chronome¬ 
ters, only measure relative distances or times, which are 
integrated quantities depending only on the dimension¬ 
less Hubble parameter H{z)/Hq, where Hq is the Hubble 
constant. Thus, provided that the equation of state of the 
total dark sector be kept constant (i.e. the sum of CDM 
and DE), distance measurements are completely degen¬ 
erate to any changes in relative composition between the 
two species. This fact was named the dark degeneracy 
in Ref. [17] and it implies that background measure¬ 
ments are incapable of measuring the cold dark matter 
(CDM) fraction Hco absent a choice of parametrization 
for wx- Measurements of large-scale structure, such as 
the CMB anisotropies, can break this degeneracy, but 
require a model for the behavior of dark-energy pertur¬ 
bations to be specified [18, 19]. One should really think 
of the density fraction of DM as a perturbation variable 
and not a background parameter [20]. While wx is fre¬ 
quently parametrized as constant and the best-fit cos¬ 
mology prefers values close to wx = —1 with a good 
determination of Oco [2], the indubitable fact is that our 
universe seems to be close to one which has a constant 
pressure at late times, making a DE component with a 
constant energy density just one of the possibilities. As 
we will show here, this is not just an academic discussion. 

Our only modification to the concordance model is to 
use a perfect-fluid dark energy with a constant (rest- 
frame) sound speed Cs and a particular designer choice of 
equation of state wx that ensures that the background 
expansion history exactly mimics ACDM (i.e. our DE has 
constant pressure). 


A. Background 


In order to maintain a fixed expansion history, we pick 
the density fraction of dark energy Q,x such that its sum 
with the CDM density fraction remains the same as the 
sum of the density fraction of the cosmological constant 
and some reference CDM density in ACDM. We 

will refit the data for the values of all the parameters in 
the paradigm of our extended model, but for the moment, 
one should think of it as the standard Planck value for 


Hx(a) + f^c(a) = f^A(a) + H“(a), (1) 


where a is the scale factor. This relation implies that the 
equation of state of the DE satisfies 


1 -I- wx(a) 


AHco 
AflcO + 


( 2 ) 


where AHc = and where the subscript 0 

denotes the value today. AHco is the only new parameter 


introduced at the level of the background and it replaces 
the more typical constant parametrization for wx- 

Generically, a violation of the null energy condition 
leads to instabilities. These are inescapable for NEC- 
violating perfect fluids (either ghost or gradient instabil¬ 
ities) [21], but may even exist in more general cases [22]. 
We thus only consider Aflc > 0, i.e. wx > — 1. 

We note here that such choices for wx as Eq. (2) 
are badly captured by the prevailing parametrizations 
in use in the community, such as a constant or the CPL 
parametrization [23, 24], which is one of the origins of our 
results. It is much more similar to models with a transi¬ 
tion in the equation of state (e.g. Ref. [25]) or Early Dark 
Energy (EDE) parametrizations (e.g. Ref. [26]), since the 
DE component tracks the CDM before the acceleration 
era. Other models with similar w(z) include quintessence 
coupled to neutrinos [27] or the generalized Chaplygin 
gas (GCG) [28]. 

The evolution of the background is essentially miss¬ 
ing all information about the microphysical content of 
the DE model, which is defined through the behavior of 
perturbations. For example, modeling the GCG pertur¬ 
bations in various ways gives rise to completely different 
constraints [29, 30]. We focus on defining the evolution 
of perturbations next. 


B. Linear properties 


We start off with a Friedmann-Robertson-Walker 
(FRW) metric with small scalar perturbations, with the 
gauge fixed to be Newtonian: 

ds^ = — (1 -I- 24')df^ -I- a^(t)(l — 2^)Sijdx^dx^ , (3) 


keeping only the scalar perturbations. 

By construction, the scalar sector of the linear pertur¬ 
bation equations for the dark energy takes the standard 
form for a perfect fluid which we do not give here (see e.g. 
[31-33]). Since the background expansion history mim¬ 
ics ACDM, the pressure of the dark energy is constant 
and therefore the adiabatic sound speed = 0 exactly. 
Therefore the only new parameter that we need to spec¬ 
ify is the sound speed Cg which we take to be constant 
and to lie in the range 0 < < 1. This is the second 

and final new parameter of our model. We can define the 
(physical) Jeans scale as 




H{z) 

(1 -I- z)Cs ’ 


( 4 ) 


which we will demonstrate separates two different 
regimes of evolution for perturbations. 

Provided that the dark energy consists of a single de¬ 
gree of freedom, and therefore has no additional internal 
freedom, the dynamics of the perturbed cosmology can 
be rewritten in a manner much more illustrative than the 
usual presentation with fluid variables (see [19, Eq. (21)]): 
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one combines the various constraints given by the Ein¬ 
stein equations to eliminate the DE fluid variables and 
to form a single equation of motion for the Newtonian 
potential which is coupled to the other matter species, 

^+3H^ + H'^+ {3H^ -h 2Er) — $ = (5) 

— Cg I 2 2 PmjHVrn I + 2 

^ = PmTTin (6) 

where the fluid variables with subscript “m" describe the 
total perturbation of the matter sector, including radia¬ 
tion, baryons, CDM and neutrinos, but not the DE. This 
system is closed by supplying the evolution equations for 
the total matter energy-momentum tensor, whether by 
assuming the matter be a collection of fluids, or by solv¬ 
ing some sort of Boltzmann hierarchy. 

Equation (5) should be contrasted with its version in 
ACDM: 

St) 

$ -f 3Ef$ + H^+ + 2H)^ = . (7) 

There are three important differences between Eqs. (5) 
and (7): 

• the presence of a sound-horizon (Jeans) term 

• the sourcing of the potential by the matter pertur¬ 
bations with a coupling c^; 

• the potential $ in our model is in fact an indepen¬ 
dent degree of freedom with its own initial condi¬ 
tions.^ In ACDM, it is completely constrained and 
not independent. 

Note that the cosmological horizon does not appear as 
a relevant scale here at all; its sound horizon is the only 
scale relevant for the DE. 

At scales larger than the DE sound horizon, k kj, 
the quasistatic approximation for the DE cannot be used 
[35]. The Jeans term on the left-hand side of Eq. (5) is 
negligible and the difference with respect to the ACDM 
case is the new source term coupled with Cg on the right- 
hand side. It affects the evolution of the potentials when¬ 
ever it is comparable to the pressure perturbation term 
Spm coming from radiation, essentially providing a new 
early integrated Sachs-Wolfe contribution. For example, 
around the time of decoupling, the CDM and radiation 
density fractions are comparable, so the new term is neg¬ 
ligible when Cg <C 1/3. For such DE sound speeds, the 
evolution of perturbations in the metric is the same as in 


^ Note that in tracker models of DE, the isocurvature modes decay 
away [34]. 


ACDM, irrespective of the value of Aflc, he. how much 
CDM has been removed and replaced with the DE. 

This means that, outside of its sound horizon, the DE 
clusters just like CDM since it cannot react causally 
to create a pressure to arrest the dustlike collapse 
along geodesics. The fact that quintessence tracking 
models behave in this way at large scales was shown in 
Ref. [34] and this is also one of the implications of the 
separate-universe approach [36]. Essentially, the part of 
the DE energy density that is tracking contributes in the 
same way as the CDM dust at these scales and therefore 
the model looks like ACDM there. 

At scales inside the Jeans scale, k ^ kj, the Jeans 
term makes the evolution equation (5) stiff, leading to 
the solution usually known as the quasistatic solution. 

When the pressure is irrelevant, this recovers the stan¬ 
dard Poisson equation of ACDM, with the gravitational 
potential driven purely by the comoving density fluctua¬ 
tions of matter. At these scales, the DE reacts causally 
and sets up a pressure profile to prevent its clustering. 
The gravitational potential is sourced only by the stan¬ 
dard species. Thus picking a AQc > 0 to replace a part of 
CDM will reduce the depth of the potentials and there¬ 
fore also the growth rate on scales inside the DE Jeans 
length. 

When the total pressure perturbation is not negligible 
compared to the total density perturbation, e.g. at de¬ 
coupling, the final term in Eq. (8) could in principle con¬ 
tribute for small DE sound speeds. On the other hand, 
we have to bear in mind that the angular scale of the 
DE Jeans scale in the CMB is proportional to the DE 
sound speed. We can use the Limber approximation to 
estimate that the multipoles at which the approximation 
(8) is valid is 

^ ^Jeans ^ ^peak 5 (9) 

Cs 

where £peak ~ 200 is the multipole at which the first 
CMB peak is centered and therefore the solution (8) is 
only possibly relevant to the observed CMB anisotropies 
for cl > 10-2. 

We illustrate the effect of our model on the CMB TT 
power-spectrum in Fig. 1, where it can be seen that 
already for Cg = 0.001, the modification of the ACDM 
curve from the terms discussed above is negligible when 
the effects from the late universe are not included. This 
can be verified by comparing a perfect Planck forecast 
with no late-time effects included for ACDM with our 
model. 

Turning to the late universe, as the mode crosses the 
DE sound horizon, the standard ACDM-like solution for 
the potential <i> valid super-Jeans is modified to (8) as 
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FIG. 1. The effect on the unlensed CMB anisotropy TT 
power spectrum of replacing roughly a third of the CDM 
(AflcO = 0.1) with our DE, as a function of its sound speed 
Cs- For clarity we show the difference between the power 
spectra for a given and the ACDM model (c^ = 0), ie 

=0). A large sound speed in¬ 
creases the amplitude of the peaks, since the DE and radia¬ 
tion Jeans horizons are close and the decay of perturbations is 
much more rapid than in ACDM. By comparing to a perfect 
Planck forecast with no lensing, we have verified that the un¬ 
lensed power-spectrum stops being sensitive to DE for sound 
speeds as large as Cg ~ 10“®. 

the DE perturbations are erased. This then modifies the 
growth function for that mode at all subsequent times, 
with the matter perturbation evolving according to 

K + 2HK-^H^il-nx)S^, = 0. ( 10 ) 

This is essentially the same equation as in the ACDM 
case, apart from the fact that even during matter domi¬ 
nation rix ~ Aflco/flco 7^ 0, which changes the growing 
mode evolution during matter domination from oc 

a to 

(5ni OC , p^l - , (11) 

5 

for fix ^ 1. This means that the matter growth func¬ 
tions at some redshift z during matter domination re¬ 
ceives a scale-dependent correction: 

G'(fc, z) cx (l-|-z)“\ fc > fcj (12) 

where z is the redshift of observation. Therefore, even 
for a fixed amount of dark energy Aflco, a larger sound 
speed means that the pivot of the suppression effect lies 
at smaller k and therefore the reduction in power is higher 
for any particular mode inside the sound horizon (see 
Fig. 2 for the effect on the CDM power spectrum or Fig. 3 
for an equivalent presentation of how the CDM power 
spectrum amplitude erg today depends on the amount 
and sound speed of our DE). 

We remind the reader that the CDM and gravitational 
potential power spectra are not affected for fc < fcj, since 
the dark energy is clustering and the potential has the 


full depth as in ACDM.^ 

In principle, the amplitude of erg is closely related to 
the rate of formation and mass function of halos. How¬ 
ever, we should note here that two different CTg’s can be 
defined: one given by the autocorrelation of the CDM 
density contrast, while the other defined from the auto¬ 
correlation of the gravitational potential. In the ACDM, 
Cg —> 0 limit, they are identical. But whenever the DE 
Jeans length lies at scales larger than 8 comoving Mpe, 
the two are different, with the former larger. On one 
hand, the forming clusters only know about the gravi¬ 
tational potential, since the CDM does not interact in 
any way but through gravity. On the other hand, the 
clusters form in configuration space and locally only the 
mass in the CDM perturbations is available to form them. 
Which of the definitions is more appropriate requires fur¬ 
ther study outside of the scope of this work and we will 
use the more usual CDM definition henceforth. 

The scale-dependent reduction of the amplitude of the 
gravitational potential affects lensing. We plot the lens¬ 
ing potential in Fig. 4, showing that the suppression is 
very strong for large sound speeds and for large Aflco- 
The effect of lensing on the CMB was detected by Planck 
by two methods: the smoothing of the power spectrum 
at small scales and using the trispectrum; and it matches 
that predicted by ACDM closely (although the power- 
spectrum method sees a 2cr excess of lensing in the con¬ 
text of ACDM) [2]. As we have already mentioned and 
show in Fig. 5, the CMB at decoupling is only mildly 
sensitive to DE sound speeds. It is the reduction of the 
CMB lensing effect in the late universe that breaks the 
degeneracy between the DE and CDM when Cg > 10“^. 

Reference [37] proposed that the integrated Sachs- 
Wolfe effect could be used to constrain the sound speed 
of dark energy and shows that WMAP data imply 
Cg < 0.04. Since the evolution of gravitational potentials 
at scales k kj in our model is purely determined by the 
evolution of the cosmological background, the integrated 
Sachs-Wolfe (ISW) effect on the CMB anisotropies also 
matches that of ACDM for a low-enough sound speed 
and this does not provide an interesting constraint. Sim¬ 
ilarly, cross-correlations of ISW with galaxy clustering 
are unlikely to constrain this model beyond some mini¬ 
mal sound speed [38]. For lower sound speeds, the above 
probes only test the expansion history.^ The novelty of 
Planck with respect to WMAP in this context is that 
CMB lensing provides new constraints which push the 
limits on Cg much further than large-scale measurements. 


^ At very late times, the comoving horizon shrinks as a result of 
the acceleration. Thus modes can in principle exit the DE sound 
horizon and have their growth increase back to ACDM rates. 

® As was argued in Ref. [36], the presence of gravitational slip 
changes the evolution of the potential at large scales. Thus the 
late-time ISW effect could still be used to constrain it, even for 
a fixed expansion history, and therefore test for modifications of 
gravity [39]. 
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FIG. 2. The effect on the CDM power spectrum at z = 0 
of replacing roughly a third of the CDM (ADcO = 0.1) with 
our dark energy. For k 3> Ho/cs, the DE does not cluster 
and causes the CDM growth rate to decrease after the mode 
crosses the DE sound horizon. This suppresses the power 
spectrum at small scales in a scale-dependent manner. Scales 
larger than the sound horizon are not affected and match 
ACDM. 



FIG. 3. The effect on the normalization of the CDM power 
spectrum, crs, at z = 0 of replacing part of the CDM with 
our dark energy. For < 10~®, the power spectrum is sig¬ 
nificantly suppressed only on scales smaller than 8 Mpc and 
therefore erg is not affected by such low sound speeds. As Cg 
increases, the power spectrum is erased on increasingly large 
scales, reducing erg. 


FIG. 4. The effect of replacing a third of the CDM (ADcO = 
0.1) on the projected gravitational lensing potential power 
spectrum as a function of the DE sound speed squared 
Cg. Large values of the sound speed erase the contribution to 
the potential from the DE even at the largest scales. As the 
sound speed is decreased, the pivot of the damping moves to 
higher multipoles, eventually giving the same lensing effect as 
ACDM. The Planck temperature power spectrum constrains 
the amplitude of the lensing signal to about 10%, which allows 
CMB lensing to break the dark degeneracy for Cg > 10“®. 

suits if they are performed at scales larger than fcj and 
will differ if they are centered on smaller scales. For 
very low DE sound speeds, there is total degeneracy. 
The scale-dependent growth rate resulting from the Jeans 
horizon allows us to alleviate some of the tension of CMB 
measurements when the sound speed is dialed up. 

It is important to stress that at no point is the effective 
Newton’s constant, describing the response of the gravi¬ 
tational potential to the CDM, smaller than its standard 
value. Deep inside the Jeans scale, they are exactly equal. 
Thus we have achieved a reduction of the growth rate 
without introducing some sort of gravity-like repulsive 
interaction which would in all likelihood be pathological. 


C. Microscopic model 


The picture we have therefore built up is that, for small 
enough DE sound speeds Cg, the dark energy perturba¬ 
tions behave as dust outside of their Jeans scale, mak¬ 
ing up for the removed CDM and giving essentially the 
same predictions as ACDM. Inside the Jeans scale, the 
DE reacts and develops a pressure which arrests its col¬ 
lapse, allowing only the CDM to cluster. This reduces the 
growth rate at small scales when compared to large, in¬ 
troducing a new scale dependence in all the observables. 
In the limit Cg —> 0, the dark energy behaves as dust 
at all scales, and therefore is completely degenerate with 
ACDM, despite its equation of state (see e.g. the Angel 
Dust model of [40] or mimetic dark matter of [41, 42]). 

As a result, the constraints that we obtain from the 
CMB are always an upper bound on the DE sound speed. 
Late-universe measurements will match the ACDM re- 


We have purposefully de-emphasized the discussion of 
the precise model of dark energy which might produce 
such phenomenology as above. Essentially any DE/MG 
model which comprises a single degree of freedom and has 
an energy-momentum tensor of perfect-fluid form gives 
this phenomenology when the solution has the equation 
of state (2) and a constant sound speed for the perturba¬ 
tions around it. Specifying the model allows for a calcu¬ 
lation of any nonlinear effects and for an estimate of the 
domain of validity of our linear model. 

As an example, the class of k-essence models [43] pro¬ 
vides for such behavior. The action is given by 

S^ = I d^xV^K{X, cj)), X = - . (13) 

The equation of state and sound speed by this dark en- 







6 


ergy are given by 

y, = _ - _ c2 = _ ^ _ (14) 

2XK^x-K^ " K^x + 2XK^xx ^ ’ 

and given some initial conditions for (/),</) a suitable K 
can be found by integrating the above. An appropriate 
choice could in fact be 

(i+cg) 

K = ¥{(/)) + M'^iX/M^) (15) 

whenever the sound speed is small Cg <C 1 and where 
V{(f)) is some sufficiently slowly varying potential, M is 
a constant mass scale. 

Another way of thinking about this model, which is 
completely equivalent at the level of this paper, is that 
we are introducing a subcomponent of the dark matter 
which exists in a condensate with a nonvanishing sound 
speed [44]. In such an interpretation, the acceleration is 
still being driven by a cosmological constant. 

One should, however, be mindful when taking such 
scalar-field superfluid models literally: the behavior of 
perturbations when nonlinear deviates from standard flu¬ 
ids, since the scalar field cannot carry vorticity and there¬ 
fore cannot virialize. The perturbations in the DE can 
become nonlinear when Cg < 10“® [45]. At the same 
time, for sound speeds Cg < 10“^ the pressure support 
that the DE can provide is never quite enough to arrest 
the collapse and such a fluid would continue to collapse 
either until the breakdown of the effective description or 
would form a black hole. 


III. RESULTS 

For our analysis we have slightly modified the 
CAMB/CosmoMC public codes [46, 47] to properly 
model the evolution of a perfect fluid with a time- 
dependent equation of state (2) and a constant sound 
speed Cg. Unless stated otherwise, we only let the param¬ 
eters AUcOj Cg and 0 mc vary, keeping the others 

fixed to their ACDM Planck 2013 best-fit values [1] . Also, 
we should mention that we scan the parameter space in 
terms of log]^o(Cs) instead of simply Cg, in order make it 
easier for the MCMC code to sample the subtle effects 
of the sound speed and recover the degeneracy for low 
Cg. For the final combined analysis, we have freed all 
parameters. 

Whenever possible, we have tried following the relevant 
parts of the analysis of [2], although we were only able 
to use the Planck likelihoods from the 2013 data release 
[1]. We first discuss the constraints from the different 
data sets in turn and then perform a combined analysis 
in Sec. IIIC. 

In order to provide a reasonable representation of the 
degeneracies, we are always including distance data to¬ 
gether with each of the perturbation-related data sets. 
Therefore, we have included the BAO measurements 


from CMASS and LOWZ of Ref. [48] , the 6DF measure¬ 
ment from Ref. [49] , the MGS measurement from Ref. [50] 
and the Union 2.1 SNe la catalog from Ref. [51], all read¬ 
ily available in the CosmoMC code. We do not include 
any measurements of the Hubble constant Hq, apart from 
a uniform prior 0.4 < h < 1.0. 

The setup most similar to ours was previously inves¬ 
tigated in [32], albeit the parametrizations for the equa¬ 
tions of state of DE used were different. At the time, 
WMAP did not provide strong constraints on the DE 
sound speed. As we show here, this is no longer the case 
in the context of a fixed background expansion history. 

A. Constraints from full CMB 

The Planck papers on cosmology [2] and on dark en¬ 
ergy and modified gravity [7] demonstrated that ACDM 
is a good fit to the data. Given the discussion in Sec. IIB, 
we expect to recover a perfect degeneracy between Uco 
and AUcO when the sound speed is sufficiently low. 

For this fit we are using both the CMB tempera¬ 
ture and low-multipole polarization power spectra from 
Planck including of course the lensing of the CMB. We 
have not included the lensing information extracted from 
the temperature trispectrum. 

As expected. Fig. 5 shows that CMB anisotropies con¬ 
strain only the sum Uco + AUcO = 0.271 ±0.004 whenever 
Cg < 10“^. For small admixtures of DE, much higher 
sound speeds are also allowed. The CMB does not show 
any preference for having this extra component, but it 
cannot rule out its existence either. 

As we have previewed already in Sec. IIB, this con¬ 
straint on Cg actually comes mainly from CMB lensing, 
which takes place in the late universe. We have gener¬ 
ated a perfect Planck forecast without lensing for ACDM 
and our model and have found that the dark degener¬ 
acy cannot be broken by the CMB power spectrum for 
Cg > 10“^. The CMB at decoupling is not very sensitive 
to clustering properties of the dark matter. 

B. Weak lensing shear 

As was noted by Planck [2, 7], the data from weak 
lensing shear are slightly incompatible with the Planck 
ACDM analysis: the preferred amplitude of the lens¬ 
ing potential found by CFHTLeNS at angular scales 
1.5' < 9 < 37.9' is lower than that expected when the ini¬ 
tial normalization of perturbations implied by the CMB 
As is extended to smaller scales and evolved forward in 
ACDM [52] . If As is kept constant at the Planck best-fit 
value, and the distance measurements are in the like¬ 
lihood, the optimal solution for WL has a lower density 
fraction UcO + AUcO = 0.244±0.084 and a reduced ampli¬ 
tude for the gravitational potential sourced by a matter 
distribution with ug = 0.761 ± 0.024. Provided that the 
sound speed of our DE Cg < 10“^, which ensures that 
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FIG. 5. The constraints on fico, AflcO and Cg from the CMB anisotropy power spectrum (including CMB lensing) and distance 
data (SNe+BAO). Provided the sound speed be cj < 10“®, the model is degenerate with ACDM for these data. Larger DE 
sound speeds are allowed if the admixture of DE is sufficiently small. The constraint on the sound speed is essentially completely 
determined by the effect of CMB lensing in the late universe on the observed CMB power spectrum, and not on the shape of 
the CMB power spectrum at decoupling. 


the data lie at scales outside of the DE Jeans horizon, 
the WL data cannot break the dark degeneracy. 

However, a new, equally good, solution for the lensing 
data appears in the presence of a small amount of our 
DE (Aflco = 0.022 ± 0.016) with a large sound speed, 
Cg > 2.5 • 10“^ which also allows for a much larger value 
of the CDM density fraction Hco = 0.265 ± 0.055. In 
such a scenario, the DE clusters only at largest scales 
and therefore the growth rate of the CDM is reduced ac¬ 
cording to Eq. (12). This reduces the amplitude of the 
lensing potential at late times despite the large flco- We 
have verified that the existence of this solution does not 
strongly depend on the nonlinear completion used for the 
matter power spectrum: just as Planck, for the fits in this 
paper we use Halofit with the parameters of Takahashi 


et al. [53]. However, we have verified that the constraints 
are not biased by comparing them with the ones ob¬ 
tained using the original Halofit parameters [54] and even 
by switching off the nonlinear correction to the matter 
power spectrum completely. We have also used the cor¬ 
relation functions from the alternative CFHTLeNS anal¬ 
ysis of Ref. [55] , verifying that it would not significantly 
influence our results. 

In principle, one could have also varied the initial fluc¬ 
tuation amplitude Ag which would have given a much 
weaker constraint for flco, compatible with the CMB fits. 
However, the CMB measures the initial amplitude very 
well and therefore WL does not significantly correct the 
posterior for it. Thus the posteriors presented in Fig. 6 
represent better the effective contribution of WL to the 













FIG. 6. The constraints on fico, AflcO and Cg from CFHTLeNS weak-lensing data and distance data (SNe+BAO). WL data 
prefer a lower amplitude of perturbations than the CMB. Since we have fixed the initial inflationary amplitude to the Planck 
A best-fit value, the solution achieves a good fit by reducing the preferred value of CDM density fraction flco compared to the 
CMB in figure 5. The dark degeneracy is only recovered for sound speeds Cg < 10~^, compared to 10~® in the case of the CMB, 
since the weak lensing data lie at much smaller scales than the CMB lensing. 

The WL data also permit an alternative, equally well-fitting, solution: flco can be larger provided that there is also a significant 
amount of DE, AflcO > 0. This DE must then have a large sound speed which ensures that its perturbation are erased already 
at large scales and the gravitational collapse of CDM proceeds more slowly at most subhorizon scales, according to the growth 
function (12). Essentially all the points in the posterior presented in the right panel which have a large sound speed come from 
this alternative solution. 


posterior of the combined data. We of course free all the 
parameters for the final combined fit. 


C. Combined constraints 

For the combined analysis, we include all the CMB 
data in the Planck anisotropy power spectrum and the 
CFHTLeNS weak lensing shear. We free all the standard 
ACDM parameters as well as AHco and Cg. 


Like in previous cases we find that a dark degeneracy 
is recovered in the combined fits for Cg ^ 10“^, i.e. in the 
parameter range as determined by the WL measurement. 
This parameter range is equivalent to the concordance 
cosmology and offers no improvement in the combined 
fit. 

However, there exists an alternative, superior solution 
already mentioned in Sec. HIB, which fits the combined 
data better than ACDM with a = 2.4 while con¬ 
taining two extra parameters. In order to discuss the 
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FIG. 7. The constraints on AflcO and cl when using the combination of Planck CMB power-spectrum data and distance data 
(BAO-|-SNe), together with weak lensing shear data from CFHTLeNS. We have imposed a hard prior cl > 10“^ in order to 
remove the infinite parameter space where our model is completely degenerate with ACDM for all the data under consideration. 
The data prefer a small admixture of DE with AflcO = 0.017 ± 0.019 and a sound speed of log]^Q(Ca) = —4.394 ± 1.779. This 
gives a cosmology which for the purposes of the CMB looks like the standard ACDM model, but the clustering is reduced at 
the scales probed by galaxy shear. 


posterior of this solution, we impose a prior to eliminate 
the part of infinite logarithmic parameter space where 
our model is completely degenerate with ACDM, cutting 
the chains at > 10“^. We show this part of the pa¬ 
rameter space in Fig. 7. 

Given this cut, the combined data prefer flco = 0.239± 
0.020 together with an admixture of Aflco = 0.017 ± 
0.019 of dark energy with a sound speed of log^g(Cg) = 
—4.394 ± 1.779. The combined density fraction flco + 
Aflco is consistent with the Planck 2013 best fit of flco = 
0.265 ± 0.011. At small scales, the DE stops to cluster 
and recovers the preferred lensing density fraction. 

We note that in this alternative solution, the amplitude 
of the CDM power spectrum is as = 0.79±0.02, which lies 


within 1(7 of the value as implied by clusters in Ref. [3], 
However, the power spectrum is quite modified at small 
scales and the predictions for nonlinear structures might 
not be related to the value of Ug in the same manner as 
in ACDM. 


IV. DISCUSSION AND CONCLUSIONS 

There is a weak but chronic tension between the CMB 
measurements of cosmological inhomogeneities and those 
obtained from the late universe, which consistently sug¬ 
gest that there is less power at late times than expected 
in ACDM. There are many systematic effects that could 
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be playing a role in biasing our interpretations of the 
low-redshift universe (baryonic physics, intrinsic align¬ 
ment, badly modeled nonlinear physics), but one should 
ask whether such a phenomenology, if it were to persist, 
could be evidence of some sort of dynamics in the dark 
energy/gravity sector. 

The usual approach to solve this is to change the geom¬ 
etry of the universe at late times (i.e. change w). How¬ 
ever, it is difficult to create a large effect on perturbations 
without significantly altering distances. The constraints 
from BAOs and SNe are now tight enough that it is very 
difficult to create a large effect consistent with these back¬ 
ground data and the CMB. Moreover, the CMB lensing 
and the WL are driven by the same physics and the ge¬ 
ometrical kernels lie close together at lower redshifts, so 
that it is difficult to produce a significant discrepancy in 
the clustering seen by these two probes with the help of 
a purely time-dependent effect. 

Taking account of the fact that the data in tension 
(CMB and WL/clusters) probe different scales, we ad¬ 
vocate an alternative approach. We keep the expansion 
history fixed and use the DE to make the CDM cluster 
differently at small and large scales. We have used a 
simple perfect fluid to model this effect. Yet this type of 
model is relatively poorly explored as a result of the fo¬ 
cus on particular classes of parametrization of w (CPL). 
We have instead exploited the dark degeneracy to freely 
vary the CDM density fraction flco while keeping H{z) 
fixed. 

The sort of parametrization we have employed gives a 
non-negligible contribution of dark energy at early times. 
This does not deform the CMB in an unacceptable man¬ 
ner, because our model strongly violates the quasistatic 
approximation typically used in the DE/MG calcula¬ 
tions. Perfect fluids at scales beyond their sound hori¬ 
zons evolve as dust, irrespective of their equation of state. 
Thus the growth rate at largest scales is completely deter¬ 
mined by H{z) and perturbation evolution in our model 
is degenerate with ACDM at large enough scales. This 
is a known and generic result for any model with no 
anisotropic stress [36] and exploited in the parametriza¬ 
tion of Ref. [56]. 

That being said, CMB lensing, detected in the Planck 
power spectrum at lOcr, puts a very strong upper limit on 
the DE sound speed. Since the CMB lensing is consistent 
with ACDM, the DE must cluster as dust on the scales 
probed by it. Then there needs to be a rapid transition in 
the clustering properties, so that the CDM growth rate 
is reduced at the smaller scales probed by WL. 

We stress that despite the fact that the CMB is a 
high-redshift observable, it constrains our model mostly 
at low redshifts. In total, the combination the com¬ 
bination of CMB and WL data fit the model better 
than ACDM with a = 2.4. This solution has a 
flco = 0.239 ± 0.020 with the presence of a small admix¬ 
ture of DE Aflco = 0.017 ± 0.019. This keeps the CMB 
anisotropies unchanged while reducing the amplitude of 
the gravitational potential inside the Jeans horizon and 


bringing the two data sets together. 

Our improved solution has reduced CTg = 0.79 ± 0.02, 
which is closer to the result reported by Planck in 2013 
from clusters (erg = 0.75 ± 0.03) [3]. This decrement in 
power continues at small scales, which could in principle 
help alleviate some of the tension between N-body sim¬ 
ulations of CDM and the observations (e.g. the missing- 
satellites problem [57] and Too Big to Fail [58]). It 
is interesting to ask whether the nonlinear behavior of 
CDM in this scale-dependent model could be mapped 
from standard ACDM N-body simulation using a method 
such as that proposed in [59] . Since in our model at small 
scales gravity is not modified and there are no new screen¬ 
ing effects, this might give a simple method for adjusting 
calibration in e.g. Halofit for this subclass of DE models 
and therefore to predict small-scale lensing and cluster 
formation with a better accuracy than we can currently. 

We have not used constraints from redshift-space dis¬ 
tortions on the CDM growth rate /erg. On one hand, we 
find that the uncertainty in the most precise measure¬ 
ments to date [60-62] is still high enough not to affect 
the posteriors significantly in the vicinity of the best-fit 
region of our scenario. More importantly, we are dealing 
with a model in which an evolving scale-dependence of 
the matter power spectrum is a key feature. Interpreting 
the growth-rate data in such a scenario is somewhat sub¬ 
tle since the growth rates are scale-dependent [63, 64] and 
we will give this question the consideration it deserves in 
a separate work. 

We should stress that we have exploited a completely 
generic feature of modified gravity. The response of the 
gravitational potential to the CDM perturbations de¬ 
pends on scale in all models apart from ACDM, whether 
because such a transition is already explicit in the 
quasistatic approximation (e.g. f{R) gravity where the 
Compton mass provides a scale [65], or see Ref. [19]) 
or because the quasistatic approximation fails beyond 
the Jeans scale (in the scenario discussed here; also see 
Ref. [35]). This implies that a scale-dependent modifica¬ 
tion of the matter power spectrum is completely generic 
in dynamical dark energy/modified gravity models and 
typically will take place in a range of modes correspond¬ 
ing to the Jeans/Compton-mass scale today and mode 
corresponding to the Jeans scale at the time that the 
DE/MG started contributing to the gravitational poten¬ 
tial significantly. 

Observations of large-scale structure are performed not 
only around a particular redshift, but inside a partic¬ 
ular range of scales probed by the survey. This scale- 
dependence can make the interpretation of measurements 
complicated, if the measurements are reported assuming 
a ACDM-like scale-independent behavior, instead of a 
rawer form closer to the observation. 

The phenomenology of such “cold dark energy” [66] 
models featuring low sound speeds as proposed here and 
exploiting a scale-dependent phenomenology, but in a 
more general context of effective field theory of dark en¬ 
ergy [20, 67-70] remains largely unexplored, but is nec- 
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essary if we are to build a full understanding of the con¬ 
straints data placed on the properties of the dark sector. 
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